Non-perturbative renormalization of the KPZ growth dynamics 



00 
(N 



I 

o 
o 



00 
(N 
(N 
O 
00 



X3 
O 

o 



X 
J3 



C. Castellano(i), M. Marsili^^^ and L. Pietronero'^'^) 
T/ie Mfcdus Salam" I.C.T.P. Strada CosUera 11, 1-34100 Trieste 
International School for Advanced Studies (SISSA), via Beirut 2-4, Trieste 1-34014 
Dipartimento di Fisica and INFM Unit, University of Rome "La Sapienza", 1-00185 Roma 

We introduce a nonperturbative renormalization approach which identifies stable fixed points 
in any dimension for the Kardar-Parisi-Zhang dynamics of rough surfaces. The usual limitations 
of real space methods to deal with anisotropic (self-afBne) scaling are overcome with an indirect 
functional renormalization. The roughness exponent a is computed for dimensions d = 1 to 8 and 
it results to be in very good agreement with the available simulations. No evidence is found for an 
upper critical dimension. We discuss how the present approach can be extended to other self-afHne 
problems. 
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Several models for the growth of rough surfaces appear 
to belong to the universahty class of the Kardar-Parisi- 
Zhang (KPZ) equation 
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Here the first term on the r.h.s. corresponds to sur- 
face tension, the second describes lateral growth and 
77 is white noise. The large scale properties of grow- 
ing surfaces are described by two exponents: the rough- 
ness exponent a, defined by the typical height fluctua- 
tion 5h ^ L°' over a linear distance L, and the dynamic 
exponent z For the KPZ equation, the exponents 
satisfy the scaling relation z -I- a = 2 as a consequence 
of Galileian invariance in the related Burgers equation 
It is then sufhcient to concentrate attention on a. In 
the usual field theory approach a perturbation theory is 
defined with respect to the non- linear term (A). The cor- 
responding renormalization group (RG) reveals |^ that 
the physics of the KPZ is related to a strong coupling 
fixed point, which is unaccessible by perturbativc meth- 
ods. Except for d = \, where, thanks to a fluctuation- 
dissipation theorem, an analytic solution is possible, the 
problem appears of strong coupling nature and not much 
can be done from this theoretical perspective. 

Extensive numerical simulations have been carried out 
for restricted solid on solid (RSOS) models [Q, for c? = 1 
to 7. These are discrete models which belong to the KPZ 
universality class |^ and are optimal for numerical simu- 
lations . These studies show a gradual decrease of the 
value of a, without evidence for an upper critical dimen- 
sion. On the other hand, several analytical approaches ||] 
suggest an upper critical dimension dc — 4: above which 
one should have a — 0. 

In strong coupling situations, real space methods per- 
form sometimes much better than continuum field the- 
ory. This is particularly true for self-organized non- 
equilibrium systems such as fractal growth sandpile 
models |^ and models with quenched disorder |^ where 
real space methods give quite accurate results. However, 



rough surfaces show anisotropic scaling properties (self- 
affinity) and, in such a case, it is difhcult to define the 
"internal degrees of freedom" for the block transforma- 
tion in real space. Here we overcome this basic problem 
by defining the renormalization transformation in an in- 
direct way that permits to integrate over the internal de- 
grees of freedom without the necessity of specifying them 
explicitly. This approach is rather general and it should 
be useful also for other self-afhne structures. 

Our starting point is to introduce a parametrization 
of the dynamics at a generic scale. Then we devise a 
RG where a scale transformation defines a flux in the 
parameter space of possible dynamics. The fixed point 
of the renormalization transformation defines the scale 
invariant dynamics and it describes the asymptotic scal- 
ing properties of the system. The key point lies in the 
choice of a correct parametrization of the dynamics which 
captures the physics of the problem but is also simple 
enough to allow for explicit calculations. The main fea- 
ture of the KPZ dynamics is lateral growth ||l|,||. This 
allows us to concentrate on just one renormalization pa- 
rameter which is the ratio between lateral and vertical 
growth rates. At the simplest level the method is fully 
analytical while for accurate results or d > 1 we resort 
to the Monte Carlo method to compute the renormal- 
ization functions. We obtain stable and attractive fixed 
points for any dimension from d = I to 8. From them we 
compute the roughness exponent a which results to be in 
excellent agreement with the available simulations Q up 
to d = 7. From the values of a and the stability of the 
fixed points we conclude that there is no upper critical 
dimension, at least up to d = 8. 

Let us consider a rough surface and define a block 
transformation for its profile as indicated in Fig. ^ The 
coarse grained description is obtained by a partition 
of the d -\- 1 dimensional space in cells of lateral size 
Lk = Lq ■ 2^ and thickness ~ L'^. Each cell of the 
partition is then declared to be full or empty accord- 
ing to some majority rule. With respect to the usual 
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RG approaches one should note that the definition of the 
coarse-graining transformation depends explicitly on the 
exponent a which, on its turn, is just what we should 
compute in the end. In addition the shape of the block 
variables changes with the scale because, since in general 
a ^ 1, the ratio hk/Lk will depend on the scale k. For 
these reasons the development of real space methods for 
self-affine problems has been extremely limited. 

With this coarse grained description of the interface in 
mind, let us now introduce the effective dynamics at the 
same level of coarse graining. The dynamics at scale k is 
defined in terms of transitions rates. At the microscopic 
scale these correspond to the possible growth processes 
occurring along the surface, as specified by the model def- 
inition (for example in Eden growth vertical and lateral 
growth have the same rates [^). At a generic scale Lk we 
generalize this situation in order to allow for proliferation 
in the RG. We set the rate of the vertical growth nom- 
inally equal to Lk and define a new parameter, \k, for 
the lateral growth. The growth rate for a generic scale k 
can then be expressed as: 

P[h{i) ^ + = ifc + Afc^ max[0,/i(j) - h{i)] (2) 

jnni 

The term Lk is the contribution of the vertical growth 
and the sum over neighbor block sites j counts the con- 
tribution of lateral growth. We assume that h{i) is ex- 
pressed in units oi hk and we neglect overhangs [ p^ . 

We can now address the following question: Given a 
description of the system at scale fc, defined by the pa- 
rameters (Lfc, hk, Afc) what will be corresponding prop- 
erties (£fc+i, hk+i, Afe+i) at the scale Lk+i = 2Lkl 

A standard answer to this question, in the real space 
RG framework, is to partition a configuration at a given 
scale into sub-configurations at a finer scale of coarse 
graining (thereby integrating over degrees of freedom). 
This method is problematic because of the self-affine na- 
ture of the problem. A different method, which we intro- 
duce here, is to consider a system of size L = Lk+2 and 
its descriptions at two different scales, say L/2 = Lk+i 
and L/A = Lk- Imposing that the computation of the 
same quantity in the two different descriptions gives the 
same result, leads to a RG equation relating the dynam- 
ics at scale fc -f 1 to that at scale k. The natural quantity 
to consider in the case of surfaces is the roughness which, 
for a system of size L, is given by 
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where h is the average height. This has the nice prop- 
erty, in the perspective of the RG calculation, that if one 
partitions the system in cells of size h 



i.e. the roughness of the system is equal to the aver- 
age of the roughness in the individual cells, plus the 
roughness of the average height of the various cells. 
Our strategy will then be to compute the roughness 
^1+2 — W^{Lk+2) over a hnear size Lk+2 first using 
the description at scale fc -I- 1 and then at scale fc. 

In order to do this we can first consider a system of 
two cells of size Lk+i and then a system of four cells at 
size fc. For simplicity we make the assumption of RSOS 
also for the coarse grained cells. This is actually an ap- 
proximation because, upon renormalization, a dynamics 
which is of the restricted type (|A/i| < 1) at the smallest 
scale can easily proliferate into an unrestricted one. How- 
ever with this restriction, the number of possible surface 
configurations is very small and an analytic treatment is 
possible. For the system of two sites we have only two 
configurations, that we classify by their height (ft.i,/i2)- 
The first one is the fiat one (1,1), while the second has 
one step (1,2). With periodic boundary conditions and 
the definition (Eq. (H)) of the dynamics at scale fc, the 
transition rates for the system of two configurations are 
Pi^2 = 2Lfc and P2^i = Lk + 2Afe. Indeed in Pi^2 only 
vertical growth is possible (in both cells) whereas in con- 
figuration 2 only one site can grow and it has two lateral 
contributions. The master equation for the probability 
Qi of configuration i is 



(5) 



and its stationary solution is (32 = 1 ~ Pi = 2/(3 -I- 2xk), 
where Xk — Xk/Lk- Using eq. (|^), we can now compute 
the fluctuations at scale fc + 1 



W'k+i 



Wl -f W\2, Xk) - Wl + Q2ixk)hll\. (6) 



Here W| is the fluctuation at scale Lk which is present in 
both conflgurations and W'^{2, Xk) is the roughness corre- 
sponding to the fluctuation of the average height at scale 
Lk- Clearly the flat conflguration does not contribute 
and the roughness of the second configuration is hi/A. 
We can relate Wl to hj. by observing that at scale fc the 
essential information we retain on the height distribution 
is only the discretization scale hk- This is equivalent to 
assuming that the distribution of the height fluctuations 
at this scale is a sum of two delta peaks 

P(iM4*(ik-|:) + l*(A. + |) (7) 



(3) It follows that Wl = hi /A and 
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Fi{xk) = 1 + 



i + 2xk 



(8) 



W'^{L) = W^{L/b) + W^{b) 



(4) 



This, at the moment, is a purely formal relation since Xk 

is unknown. Clearly eq. (^) also relates 

Xk+i, and combining these two relations we get: 
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Fi{xk)Fi{xk+i 



(9) 



The roughness M^|_|_2 can also be computed directly us- 
ing the dynamics at scale k. For this we need to consider 
a system of four sites, at scale k. There are, always within 
the RSOS approximation and with periodic boundary 
conditions, six possible configurations of (hi, h2, h^, hi). 
The transition rates Pi~tj are computed using eq. (||) and 
the stationary solution gi of the master equation ^ is 
easily found. This allows one to compute the roughness 
f*-"" 4 cells system following exactly the same 
procedure used for the 2 cells system, obtaining 



F2{xk) - 1 + 



51 + 86a;fc 



404 



47 + IQQxk 



68x1 
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Consistency of the two descriptions requires that Eqs. 

and ( |To| ) must yield the same quantity. This leads to 
our basic RG transformation for Xk+i as a function of Xk 



F2{xk)/Fi{xk) 



(11) 



It can also be cast into the more conventional form 
Xk+i = Rixk). A fixed point x* = R{x*) ~ 2.08779... 
exists and we can compute the roughness exponent as 



a 



W, 



k+l 



W'k 



= -log2Fi(x*) 



0.177352... (12) 



The fixed point is also attractive. Close to the fixed point, 
we have Xk+i — x* = R{xk) — x* ~ R'{x*){xk — a;*). If 
\R'{x*)\ < 1, the distance Xk — x* ^ as fc oo. We 
find R'{x*) — —0.03548... i.e. x* is strongly attractive. 

The value of the exponent is clearly too small with re- 
spect to the exact one {a = 1/2). This is due to the 
RSOS approximation which, on the other hand, allows 
for a fully analytical treatment. The point is that even 
if at the microscopic scale the dynamics of the surface is 
of restricted type (RSOS), the effective dynamics for the 
block variables at a generic scale k defined by the renor- 
malization procedure can proliferate in a non restricted 
dynamics. Allowing larger steps makes the problem ana- 
lytically intractable. Still the calculation of the functions 
Fi and F2 can be done using a Monte Carlo method: For 
each value of the parameter Xk — Xk/Lk, the growth pro- 
cess in cells of sizes b — 2 and 4 described by eq. (|^), is 
numerically simulated. -^1(2;) and -^2(2^) are then com- 
puted as time averages in the stationary state from the 
roughness W^{b). In this way one sees that, allowing for 
larger steps, the value of a approaches the exact value 
a = 0.5 with great accuracy. Convergence is reached 
when steps \Ah\ « 8 are allowed. Going further, i.e. to 



16, a does not change any more |11|. 

The Monte Carlo calculation of the functions Fi{x) 
and ^2(2;) also allows to deal with higher dimensions 
d > 1. The function Fi is then related to the roughness 



of a cell of size 2"^, while F2 refers to 4''. The method is 
exactly the same as for the one dimensional case. The 
results for d ~ 1, . . .8 are reported in Fig. Instead 
of R{x), Fig. ^ displays the function R{a) where a(x) is 
given in the second equality of eq. (|l2|). This method 
is graphically convenient since in this way the exponents 
can be read from the intersection of R{a) with the line a 
and the stability of the fixed point can be inferred from 
the slope R'{a*) of the function at the fixed point. We 
see that for d = 1, . . . , 8 a fixed point with a > exists 
and it is attractive (see inset). The values of the expo- 
nent a are reported in table , together with the fixed 
point value x*{d) for d = 1, . . . , 8. For comparison, table 
also reports the exponent a obtained in ref. |Q] from nu- 
merical simulations up to d = 7. The agreement between 
our RG results and the simulations is rather impressive 
and it provides a strong support to the validity of both 
methods. 

The extrapolation to d ^ 00 of R'{x*) suggests that 
the fixed point never becomes unstable, leading to the 
conclusion that no finite upper critical dimension exists 
for the KPZ dynamics (note that a = in our scheme 
would imply a stable fixed point at x* — 00). On the ba- 
sis of these results, the arguments leading to an upper 
critical dimension dc = 4 should probably be reconsid- 
ered. Furthermore our results suggests a behavior of the 
type a{d) ^ 1/d and x*{d) ^ 2'^ for large d. Indeed an 
expansion of our RG method for large dimensions, to be 
reported elsewhere ||ll[], leads to such a behavior. Finally 
we verified that our RG scheme is stable with respect to 
the change of the cell sizes used. Indeed a RG function 
Xk+\ — R(xk) can also be defined considering the func- 
tions F,+i(x) = W-feVi/I^I, F,{x) = Wl^.lWl and 
Fi(a;) (the case discussed above being £ = 1). In d = 1 
we found essentially the same results, whereas in larger 
dimensions we found the same qualitati ve p icture with a 
weak dependence of the exponent on I 111 . 

In summary we have introduced a new method to 
renormalize self-affine dynamics in real space. The key 
technical element is to effectively integrate over inter- 
nal degrees of freedom in an indirect way via a consis- 
tency relation. This allows us to overcome the difficulty 
of defining a geometric fine graining procedure for self- 
affine structures. In some sense, the method general- 
izes the Fixed Scale Transformation Q idea to self-affine 
structures. There are two key steps: first we introduce a 
description of the system and the dynamics at a generic 
scale Lk in terms of the parameters hk and = XkLk. 
Secondly we derive the RG equation Xk+i — R{xk) with 
a new indirect method. The attractive fixed point x* 
(self-organization) of the RG transformation identifies 
the scale-invariant dynamics and allows us to compute 
the exponent a. 

From the point of view of the approximations involved 
in this method, we observe that the first step is essen- 
tially a description tool and, as such, it is in principle 
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exact. The approximations arise in the identification of 
the scale invariant dynamics, which can be made more 
and more accurate with the introduction of additional 
proliferation parameters. This explains the striking ac- 
curacy of the present method as compared to usual real 
space RG methods for equilibrium phenomena. 

Usually real space RG methods are at their best for low 
dimensional systems. The topological complexity of in- 
teractions indeed increases with dimension thus requiring 
more and more drastic approximations for larger d. This 
is not the case here: Our method is exactly the same 
in any dimension and the computational complexity is 
overcome with Monte Carlo tecniques. 

The nature of our approach is quite different from 
the usual perturbation techniques for KPZ. Indeed our 
parametrization of the dynamics does not include the 
Edwards- Wilkinson equation (A = in eq. (1)). Rather 
the non-interacting limit of our theory is random depo- 
sition (A = = 0). 

Finally we observe that this general method can be 
applied to a wide variety of other dynamic critical phe- 
nomena. In this sense the present application to the KPZ 
problem is a very promising example. 

We acknowledge interesting discussions with A. Mari- 
tan, A. Stella, C. Tebaldi and A. Vespignani. 
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FIG. 1. Coarse grained description at scale of a growing 
interface. 
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FIG. 2. Renormalization group transformation R{a) for 
d = 1,...,8. The intersections with the line a yields the 
fixed point value a(x*) of the exponent. The exponents a 
and the slope R (a) at the fixed point, are plotted in the in- 
set as a function of d. Since |.R'| < 1 we conclude that the 
fixed points are stable. 
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TABLE I. Results of the RG approach compared with numerical simulations. The first line gives the value of the fixed point 
critical parameter x* . The fixed point is attractive and the corresponding exponent qrg is listed in the second line. The third 
line line reports the numerical results of ref. [4]. The excellent agreement agreement between our RG result and numerical 
simulations provides a strong support for the validity of both, and it points to the absence an upper critical dimension for this 
problem. 
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